Weak and strong wave turbulence spectra for elastic thin plate 
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Variety of statistically steady energy spectra in elastic wave turbulence have been reported in 
numerical simulations, experiments, and theoretical studies. Focusing on the energy levels of the 
system, we have performed direct numerical simulations according to the Foppl-von Karman equa- 
tion, and successfully reproduced the variability of the energy spectra by changing the magnitude 
of external force systematically. When the total energies in wave fields are small, the energy spectra 
are close to a statistically steady solution of the kinetic equation in the weak turbulence theory. On 
the other hand, in large-energy wave fields, another self-similar spectrum is found. Coexistence of 
the weakly nonlinear spectrum in large wavenumbers and the strongly nonlinear spectrum in small 
wavenumbers are also found in moderate energy wave fields. 



PACS numbers: 62.30.+d, 05.45.-a, 46.40.-f 

There exist two types of "turbulence theories" started 
by Kolmogorov and Zakharov. The former, based on the 
self-similarity and the dimensional analysis, is applied on 
vortical flows governed by the Navier-Stokes equation. 
The latter, based on the random phase approximation, is 
applied on wave fields, and is called the weak turbulence 
theory (WTT) [1]. While Kolmogorov's turbulence the- 
ory is applicable to strongly nonlinear turbulence, WTT, 
which is mathematically sophisticated, is applicable to 
weakly nonlinear wave fields @, Q • 

Elastic waves propagating in thin elastic plates are 
considered a testing ground for numerical or experimen- 
tal verification of the applicability of WTT. The ki- 
netic equation obtained in WTT has two statistically 
steady solutions One is in a non-equilibrium state, 
where the azimuthally- integrated energy spectrum £{k) 
as a function of wavenumbers k is given as £{k) oc 

1 /3 

A;(log(fc*/fc)) /0 . The other is in a thermal equilibrium 
state, where £{k) oc k. The function form of the non- 
equilibrium solution is determined so that the energy flux 
is constant, but the large-wavenumber cutoff fc* is inde- 
terminate in the framework of WTT. Additionally, WTT 
cannot determine the direction of the energy flux in the 
non-equilibrium state in the elastic wave turbulence [![. 
In these respects, the elastic wave turbulence is atypical 
as a weak turbulence system. 

During et al. obtained the energy spectrum corre- 

1 /3 

sponding to £{k) oc k (log(fc*/fc)) ' by a direct numeri- 
cal simulation (DNS). On the other hand, energy spectra 
close to £(k) oc fc~ a2 arc reported in experiments using 
thin elastic steel plates 0, Q . The difference in the en- 
ergy spectra is explained in terms of the anisotropy of the 
system in Ref. [7j, and the discreteness of the numerical 
simulation in Ref. Q. Moreover, dimensional analysis 
based on the self-similarity predicts energy spectra cor- 
responding to the energy cascade, £{k) oc fc _1 , and to 



the wave action cascade, £{k) oc A: -1 / 3 [§]. 

We optimistically believe that there should exist a sim- 
ple unified explanation for the variability of the energy 
spectra including the theoretically-predicted spectra. In 
this work, a series of DNS is performed according to a 
basic equation. Non-equilibrium steady wave turbulent 
states are obtained by adding external forces to small 
wavenumbers and dissipation to large wavenumbers. The 
numerical results provide a unified perspective on the 
variability of the spectra. 

The governing equation for the lateral displacement Q 
and the momentum p in a thin elastic plate is called the 
Foppl-von Karman (FvK) equation [lOj, 



dtp = 
A 2 x = 
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where x is t ne Airy stress potential. The Laplace oper- 
ator and the Monge-Ampere operator are expressed as 
A and {f,g} = d xx fd yy g + d yy fd xx g - 2d xy fd xy g, re- 
spectively. The Young's modulus E, the Poisson ratio a, 
and the density p are the physical properties of the plate. 
The thickness of the plate is expressed by h. The FvK 
equation ([TJ holds even for large displacements when the 
gradient of the displacement is smaller than unity. 

The linear dispersion relation between a wavenumber 
vector k and the corresponding frequency Wfc is given as 



Eh 2 



12(l-CT 2 )/7 



k 2 



(2) 



where k = \k\. The complex amplitude = {puJkCk + 
ipk) I \/2puJk represents an elementary wave, where pk 
and £fc are the Fourier components of p and respec- 
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tively. The governing equation can be rewritten as 



(3) 



where Afk symbolically expresses the four-wave nonlin- 
ear interaction terms. Note that no three-wave interac- 
tions exist in this system. (See Ref. [3] for the nonlinear 
terms.) This equation can also be written in Hamiltonian 
form, ihk = SJi/Sa^, where 5/8a* k denotes the functional 
derivative with respect to the complex conjugate of a^. 

WTT can be applied to the elastic wave turbulence, 
if the nonlinear interactions among a k are weak and 
the random phase approximation is valid. The wave 
action ny. is defined as (aka* k ,) = nkSkk', where (• • ■ ) 
expresses the ensemble average, and Skk' is Kronecker's 
delta. When the system is statistically isotropic, we ob- 
tain the relation £{k) oc kujk n k- WTT considers the 
energy transfer due to the weakly nonlinear resonant in- 
teractions among four wavenumbers. We have confirmed 
that this system has huge number of the quartets which 
exactly satisfy the four-wave resonant conditions. Then, 
the system is little affected by the numerical discretiza- 
tion and it never goes into the frozen turbulence fill ] - 

To investigate the statistically steady states, the fol- 
lowing forcing and dissipative mechanisms are added to 
Eq. ([3]) as usually done in turbulence simulations. Energy 
is inputted at small wavenumbers and dissipated at large 
ones. In the forced wavenumbers |fe| < 87r, the absolute 
value of the complex amplitude |<2fc| is kept constant as C 
by multiplying the factor CJ \ | at each time step to con- 
trol the energy of the system easily. In spite of fixing the 
absolute value, the phase of each mode evolves accord- 
ing to Eq. ([3]). The artificial eighth-order hyper- viscosity 
— zy(|fc|//cd) 8 Ofc is added to the right-hand side of Eq. (0}, 
where v is a dissipation coefficient and fed is a dissipa- 
tion wavenumber. The periodic boundary condition is 
adopted to match the assumptions in turbulence theory. 
The pseudo-spectral method with the aliasing removal by 
the 4/2 law is employed to obtain the nonlinear terms. 
The fourth-order Runge-Kutta method is employed for 
the time integration, and the linear dispersive term and 
dissipative term are implicitly solved to improve numer- 
ical stability. Independently of the initial conditions, the 
system is attracted to a statistically steady state where 
the external force and the dissipation balance. 

The material properties in the experiment with steel by 
Boudaoud et al. [5j are adopted in the present numerical 
simulations. Namely, p — 7.8 x 10 3 kg/m 3 , E — 2.0 x 
10 n Pa, a = 0.30, and h = 5.0 x 10" 4 m, respectively. The 
plate is supposed to have the periodic boundary lmx lm. 

A series of DNS with 512 2 (256 2 alias-free) modes [l2| 
is performed in nine energy levels E1-E9. The energy 
levels are set by doubling the constant C in the forcing 
term. The ratio of the stretching energy, which is the 
nonlinear part of the Hamiltonian, to the bending energy, 
which is the linear part of the Hamiltonian, represents 
the level of the nonlinearity that the entire wave field has. 




FIG. 1. Variability of the energy spectra due to energy levels. 
The energy levels are numbered from bottom to top. The 
dashed-dotted line and the broken curve respectively show k 
and k(\og(k*/k)) 1/3 where fc* = 144-k. 



The ratios are averagely 1/600 in El, 1/20 in E5, and 1/4 
in E9. Note that the FvK equation is justified since the 
root mean squares of the gradients of the displacements 
(|VC| 2 ) 1/2 are respectively 5.0 x 10~ 3 in El, 4.3 x 10~ 2 
in E5, and 4.8 x 10 _1 in E9. The standard deviations of 
the lateral displacements, (C 2 ) 1 ^ 2 , are 1.9 x 10~ 4 m in El, 
1.4 x 10~ 3 m in E5, and 1.4 x 10~ 2 m in E9. The energy 
spectra for E1-E9 are drawn in Fig. [TJ These spectra 
are time-averaged during the statistically steady states, 
and the fluctuations are too small to affect the power- 
law exponents. Note that all the spectra in |fc| < 87T are 
proportional to fc 3 owing to the forcing term. 

We observe the spectral similarity in the low energy 
levels, i.e., in E1-E3. In the small energies, that is, when 
the nonlinearity is weak, the energy spectra approach the 
spectrum predicted by WTT, £ (k) oc fc^log^/fc)) 1 / 3 , 
which fits the spectra by choosing fed = 1447T as fc*. 
As the energy levels get elevated, i.e., in E4-E6, the 
downward-sloping spectrum spreads from the smaller- 
wavenumber range. We observe another spectral simi- 
larity in the high energy levels, i.e., in E7-E9, where the 
nonlinearity is relatively strong. The energy spectra are 
estimated as £{k) oc /j-°-30±o.04 me thod of least 

squares. From now, for simplicity we use the expression 
"strongly nonlinear" for the relatively strongly nonlin- 
ear regime where WTT is not justified. For a wave field 
with the energy between E6 and E7, one may observe 
£{k) oc fc -0 - 2 as in Refs. [j| and [§[, if one forcibly fit 
the spectra to a single self-similar spectrum. Though 
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FIG. 2. The absolute values of the complex amplitudes \au\ 
during three periods, (a): Large energy E9, (b): moderate 
energy E5, (c): small energy El. 
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FIG. 3. The frequency spectrum |afe(r2)| 2 of each wavenum- 
ber. (a): Large energy E9, (b): moderate energy E5, (c): 
small energy El, and the spectrum around Q/oJk = —1 is en- 
larged in the inset. Each spectrum is normalized for visibility 
so that its maximum is unity. 



their experimental energy spectra appear to collapse to 
a similar function form, the change of the energies are 
as small as one or two decades. In the experiment on 
decaying turbulence where the change of the energies are 
relatively large Q , the energy spectra approach the pre- 
diction of WTT as the energy decays. Our simulations 
reproduce all these experimental results by changing the 
magnitude of the external force as shown in Fig. [TJ Note 
that the reproducibility supports the applicability of the 
FvK equation in all these energy levels in addition to the 
spectral similarities mentioned above. 

To investigate the nonlinearities in each scales, the 
time evolutions of the magnitudes of elementary waves, 
|dfc|, which would remain constant in the absence of the 
nonlinear interactions, are shown in Fig. [5J We select 
k = (16tt,0), (32tt,0), (64tt,0) and (128tt,0) as the rep- 
resentatives in each scales. Each evolution is drawn 
for three periods determined by the linear dispersion, 
3Tfc = 67r/cjfc. As the energy becomes larger from El 
to E9, the fluctuations of |afe|'s become larger and faster 
as the overall trend. In Fig.f^a), i.e., in the large energy 
E9, all |cifc|'s are far from constant. The nonlinear inter- 
actions are active in all the scales. On the other hand, in 
Fig.lSJc), i.e., in the small energy El, all |ofc|'s are almost 
constant in time. The nonlinearity at each wavenum- 
ber is uniformly weak. Then, WTT is applicable for El, 
but not for E9. In Fig. H{b), i.e., in the moderate en- 



ergy E5, the smaller- wavenumber waves have larger and 
faster fluctuations. It suggests that the failure of WTT, 
which appears at the wavenumbers where the linear and 
nonlinear time scales are comparable 13j, arises in the 
smaller-wavenumber range. 

In Fig. [31 we examine the frequency spectrum |afc(f2)| 2 
of <Zfc(i) at each wavenumber, where f2 is the angular 
frequency. The frequency spectra give the conjugate 
properties with the time evolution of |afc(i)| shown in 
Fig. [2] Note that if the nonlinear terms were neglected 
in Eq. ([3]), the frequency spectra would be the line seg- 
ment at O/wfe = — 1. As shown in Fig. OJa) and (b), for 
the higher energy levels and for the smaller wavenumbers, 
the frequency spectra become broader, and the peaks of 
the frequency spectrum distribution shift more to the 
(negatively) larger frequencies from the linear dispersion 
relation 14|. For the strongly nonlinear waves, the lin- 
ear dispersion relation does not characterize the system. 
When the energy is small, in Fig. [3]Jc), the peak fre- 
quencies almost coincide with those given by the linear 
dispersion relation. Even in this case the frequency shift 
to the larger frequencies in the small wavenumbers are 
found in the inset. It is consistent with the observations 
in the experiments [l5j , though the boundary conditions 
may slightly affect the deviation as they described. 
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The mean frequency shift is estimated by self interac- 
tions, which act as the dispersive term |3j. Then, Eq. ([3]) 
can be rewritten as a k = — i(uJk+Auj k )a k +Af k , where Afl. 
expresses the non-self nonlinear interactions. The mean 
frequency shift due to the self interactions is given as 



Auj k = 



E 



\k x k' 



-f 4p 2 uj k Lu k i \k — k' 



(K'| 2 + |a-fc<| 2 ) 



Suppose that the energy spectrum is self-similar as 
£ (k) oc k a and therefore |afc| 2 oc k a ~ 3 . Then, the sum- 
mation converges if —1 < a < 3, and the mean fre- 
quency shift is proportional to k a ~ x . Both for the weak 
turbulence spectrum (a = 1) and for the strongly non- 
linear spectrum (a = —0.3 shown in Fig. [1]) the mean 
frequency shift is relatively large in the small wavenum- 
bers compared with the linear dispersion, ui k oc k 2 . The 
large mean frequency shift in the small wavenumbers is 
consistent with the appearance of the strongly nonlinear 
spectra in the small wavenumbers in Fig. [TJ 

We also examine the isotropy of the system by divid- 
ing the wavenumber plane azimuthally into eight regions. 
The maximal difference between the energy in one of the 
regions and the mean energy in the eight regions is used 
as a measure of isotropy. Even in the highest energy level 
E9, the instantaneous relative difference is smaller than 
10%, and the time- or ensemble-averaging would make 
the difference smaller. Therefore, the wave field in E9 
is statistically isotropic. The isotropy is validated also 
by the two-dimensional energy spectrum. The isotropy 
in other energy levels is similarly confirmed. While the 
deviation of the power-law exponents observed in the ear- 
lier studies from the prediction of WTT is accounted for 
by the anisotropy in Ref. Q, the isotropy is maintained 
in all of our simulations. We therefore give an alternative 
mechanism for the deviation, which is accounted for by 
the nonlinearity. 

In this Letter, focusing on energy levels, we have shown 
that the level of the nonlinearity provides the unified 
perspective on the variability of the spectra in the ear- 
lier studies 0-0] • The energy spectra in the low en- 
ergy levels agree with the weak turbulence spectrum 

1/3 

£(k) oc k (log(fc„»/fc)) 1 . Less dissipative plates must be 
more weakly forced to reproduce the weak turbulence 
spectrum in experiments. The simulations also show an- 
other self-similar spectrum whose power-law exponent 
is approximately —0.30 in the high energy levels. The 
power- law exponent reminds us of —1/3 for the "inverse 
wave action cascade spectrum [la" , which follows from 
the dimensional analysis, in Ref. |9J. Note that the phys- 
ical picture of the wave action is unclear in strongly non- 
linear wave turbulence and the wave action is not con- 
served even under the kinetic equation in WTT of the 
present system owing to the 1 — > 3 and 3 — > 1 asymmet- 
ric resonant interactions. Therefore, it is still an open 
question how the energy spectra in the high energy lev- 



els are created. 

Moreover, in the moderate energy levels, the coexis- 
tence of the strongly nonlinear spectrum in the small 
wavenumbers and the weakly nonlinear spectrum in the 
large wavenumbers is found. The coexistence of the 
weakly and strongly nonlinear turbulence is predicted in 
several anisotropic wave turbulence systems [17| . In addi- 
tion, the energy equipartition in the small wavenumbers 
and the weak turbulence spectrum in the large wavenum- 
bers are simultaneously observed in a one-dimensional 
mathematical model of wave turbulence [48]. Since the 
FvK equation describes real physical dynamics, it is our 
future work to clarify the relation between real-space 
structures and the fluxes of conserved quantities. 
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